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Abstract 

The Lagrangian density of relativistic mean-field (RMF) theory with non-linear derivative 
(NLD) interactions is applied to isospin asymmetric nuclear matter. We study the symme- 
try energy and the density and energy dependences of nucleon selfenergies. At high baryon 
densities a soft symmetry energy is obtained. The energy dependence of the isovector 
selfenergy suppresses the Lane-type optical potential with increasing energy and predicts 
a p-meson induced mass splitting between protons and neutrons in isospin asymmetric 
matter. 
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1. Introduction 



The study of isospin asymmetric nuclear matter is closely related to a few issues which 
are presently the subjects of ongoing theoretical and experimental researches. At first, 
the isospin dependence of the nuclear equation of state (EoS) shows up in the density 
dependence of the symmetry energy. A quantity which drives the structure of neutron rich 
nuclei and at supra-high densities the physics of supernovae explosions and the structure 
of compact neutron stars [l|, @|. For instance, the stiffness [ij of the symmetry energy is 
directly related to the proton-fraction in /3-equilibrated neutron stars, which is limited by 
the onset of the URCA-process jij. However, the symmetry energy is empirically known 
up to saturation densities, where it can be determined in nuclear structure studies from 
Giant and Pigmy resonances and from neutron skin thickness jsj] . In the laboratory the 
symmetry energy can be explored in heavy ion collisions, where the knowledge of the 
hig h density EoS is mandatory 0, 0] for the description of the collision dynamics (si, 0, 
10 . 11 . 12 . 13 ] within transport theory 14, 15|. Different heavy ion observables, e.g., 

16] of particles and/or 
17 1 are sensitive to the 



collective transverse (in-plane) and elliptic (out-of-plane) flows 
light fragments, as well as subthreshold strangeness production 
stiffness of the EoS. Furthermore, the isospin tracing 
flows 



20, 21, 22, 23 



lif . isoscaling 




and isospin ratios of produced mesons 
plausible tools to constrain the isovector EoS at high densities and energies, 
heavy ion collisions allow to probe not onl y th e density of created matter, but also the 



collective isospin 
represent 
Moreover, 



energy dependence of the nuclear potential [28|, |29|, [30 . 
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Theoretically, the density and energy dependence of the isovector EoS beyond the 
saturation density is still controversial. In non-relativistic Skyrme models both possibilities 
of a soft and a stiff symmetry energies [H, 0, [2ll 22] can be obtained. In the RMF models 
the symmetry energy is determined by the Lorentz-structure of the isovector potential [H, 
which is described by the vector-isovector p-field, or by p- and an attractive scalar 5- 



fields [3l|, l32j. In both cases, the symmetry energy rises with increasing baryon density 



leading to a stiff isovector EoS at high densities. The models with density dependent (DD) 



couplings [33|, |34|, 1351 



Hartree-Fock (DBHF) calculations [38 



361 ] predict a soft symmetry energy [37 

'3 



while the Dirac-Brueckner- 
and other RMF models with complementary 
exchange (Fock) terms (density- dependent relativistic Hartree-Fock (RHF) [40|, |4l|) favour 
again a stiff one. 

Concerning the energy dependence, the Schrodinger equivalent optical potential in sym- 



metric matter is well known phenomenologically [42J. However, different empirical analyses 



give different predictions for the energy dependence of isovector Lane-potential [43j, |44|, |45 



46l . 147] . The microscopic DBHF [38] and BHF [48|, |49( models predict a rather constant 
optical potential at low momenta with a strong decrease at high energies. In Skyrme ap- 
proaches both cases of a strong increase and a strong decrease of the Lane-potential as 
a function of energy can be realized [2J. The DD models do not contain any explicit en- 
ergy dependence, except if one modifies the original DD formalism 33|, |34j with additional 
interaction terms 



37 



however, with the cost of additional parameters. Note that, the 
energy dependence of the optical potential has been so far a crucial problem in standard 
RMF, where the energy dependence appears only due to relativistic effects. 

Another question concerns the quasi-particle properties of nucleons in the nuclear mat- 
ter. The isospin-dependent forces may result in a splitting between the in-medium masses 
of protons and neutrons which can strongly influence the dynamics in heavy-ion colli- 



sions 



21, 22 



However, it is still under debate, whether m* 
tron rich matter. In non-relativistic approaches the in-medium particle mass is described 



< m p or m* 



> m* holds in neu- 



through the energy dependence of the real part of the phenomenological mean-field [50 
Again, the Skyrme-like interactions predict both possibilities for the non-relativistic 
in-medium mass 50] . In relativistic models the in-medium dependence of the particle 
mass arises from the Lorentz-scalar character of the relativistic mean-field, i.e., from the 
scalar selfenergy 50j. It is often denoted as "effective mass" or "Dirac mass", and these 



definitions we will use in this work. As in the non-relativistic case, same situation holds for 
the DBHF models where some calculations predict for the Dirac mass m* < m* [38J and 



others m* > m* 



51, 52 



In RMF models the mass splitting is merely fixed by the Lorentz 
structure of the isovector mean-field. Due to the Lorentz-vector character of the pNN- 
interaction the minimal crcup-models do not generate any in-medium mass splitting, except 
if one in addition introduces the Lorentz-scalar 5-field in the isovector channel 
or if one includes explicitly exchange (Fock) terms in the relativistic formalism 



31 
40 



41(. 



Including the 5-field one gets m* < m*, while in the RHF models both cases m* ^ 



are possible. So far these were the only 
relativistic mean field theory 3l|, |32|, |37|, |4 



ossibilities to generate the mass splitting in the 



41 



In this work we study the isovector interactions in asymmetric nuclear matter. Our 



2 



studies are based on the non-linear derivative (NLD) model to RMF proposed in [53j. This 
approach describes remarkably well the bulk properties of nuclear matter and, in particular, 
the density dependence of the nuclear EoS and the energy dependence of the proton- 
nucleus [53] as well as the antiproton-nucleus Schrddinger equivalent optical potentials [54 . 
The present paper is organized as follows. In Section 2 we briefly outline the NLD model. 
In Section 3 the approach is applied to isospin symmetric and isospin asymmetric nuclear 
matter. The numerical realization and details of the calculations are described in Section 4. 
The results and discussions are presented in Section 5. We finally end up with the summary 
in Section 6. 



2. Non-Linear Derivative (NLD) Model 



The NLD approach is based on the Lagrangian density of relativistic hadrodynamics 
(RHD) (55I . 56 ] and describes the in- medium interaction of nucleons through the exchange 
of auxiliary scalar-isoscalar a, vector-isoscalar uj^ and vector-isovector p M meson fields. 
The Lagrangian density is given by 



1 22 1, 

2 2 M 



F 1 " + \m 2 p P»P» 



—G fj, u G 



fJtU 



+ C 



int 



where the first two terms describe the (symmetrized) Lagrangian for the nucleon field 
\]/ = f^„, ^n) T with bare mass m and the second line contains the standard Lagrangian 



densities for a, w' 1 and p' 



fields with masses m a , m u 



and m„, respectively. The strength 



tensors are defined as F^" = d^u v — d u tu^ and G = d^p 11 — d u p fl for the vector-isoscalar 
and vector-isovector fields, correspondingly. 

In Eq. ([Q) the last term Ci n t r contains the interaction between the meson and nucleon 



fields. In the conventional RHD 55|, |56] the meson fields couple to the nucleons via the cor- 
responding Lorentz structures \I/\I/o", -^^^uj^ and -\I'7 At r\I'p A j. Such minimal Walecka-type 
interactions describe successfully the static properties of nuclear matter around saturation 
density. However, they are not sufficient for the description of dynamical situations which 
occur in nucleon-nucleus and heavy-ion collisions, where this minimal model results in a 
wrong energy dependence of the mean-field. This has been a crucial problem of the RMF 
models which already attracted much attention in the past 3(J 37, 57 1. A possible solution 
to this problem has been suggested in 28j, where momentum- dependent form factors which 
suppress the high momentum components of fields were introduced. This idea has been 
followed in the NLD model 53j where it was generalized in a manifestly covariant way. 

The NLD interaction Lagrangian introduces the non-linear derivative terms into the 
conventional RHD meson-nucleon vertices. In particular, the Lagrangian density depends 
on higher-order partial derivatives of the nucleon field ^ 



C int = C(V, d a q, d a dfl, 



tf, d a y, d a dp*, •••) 



(2) 



3 



with the symmetrized meson-nucleon interactions given by (r denotes the isospin operator) 



int 



9a 
2 

2 
9p 
2 



(3) 



with obvious notations for the fields and their couplings. The new operator T> acts on the 
nucleon fields and is the generic non-linear function of partial derivatives. In the model of 
Ref. 53| a simple exponential cut-off function has been assumed 



V-= exp 



V-— exp 




oo 

E 

n=0 



E 

n=0 



(id a v a /A)\ m/A 



rv. 



(4) 



(5) 



In Eqs. fl4]) and ([5]) t> a denotes an auxiliary dimensionless 4-vector and A stands for the 
cut-off parameter. A way to realize the meaning of the cut-off function is the presence 
of the momentum (or energy) dependent form factors in the one-boson-exchange models, 
which cut the high energy behaviour of the meson-exchange vertices 58|, |59j. Furthermore, 
in the limiting case of A — > oo the standard RHD (or Walecka) model is retained. 

The Euler-Lagrange equations and the Noether currents are derived from variational 
principles which are applied to the NLD Lagrangian containing higher order derivatives of 
the spinor fields \1/ and \1>. The generalized Euler-Lagrange equations read (</> = \?) 



TT7 - O ai 



dC 



d(d ai (j)) 



+ 9 ai d, 



dC 



a 2 



d(d ai d, 



a: 2 1 



H h 

{-) n d ai d Q 



(6) 



dC 



d(d ai d a2 ■ ■ ■ d an c 



The invariance of the NLD Lagrangian under global phase transformations results in 



the following Noether current J M |53 
dC „ dC 



J" 



-d. 



+ 

+ 
+ 



d(d n d P (, 
dC 



dC 



a 



d(dndpd y (, 

dC 
d(d n dpd^ 

dC 
d{d ll dpd 1 (i 



± 



d,s 



dC 



d{d IJ ,dpd 1 ds( 1 



± 



(7) 
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Furthermore, the invariance under space-time translations gives the energy-momentum 
tensor jSj 



dC 



dr. 



dC 



+ 

+ 
+ 



d{d^dfi(t)) 

dC 
d(dpd P (j>) 



dC 



a 



d{d ll dpd^ 

dC 
d{d lM d p d 1 i 

dC 
didpdpdtf 



± • ■ ■ 



a 



dC 



d(dpdidpde 



± 



d v (j> 



The Euler-Lagrange equations fl6]) and the Noether theorems (JTj) and OH]) contain now an 
infinite series of higher order derivatives of the nucleon field. 

It was shown in (Hi} that in mean-field approximation (up to derivatives of the meson 
fields which in any case vanish in nuclear matter) these infinite series can be resumed 
exactly. Indeed, the application of the generalized Euler-Lagrange equations ([6]) to the full 



Lagrangian density (J3J) with respect to the field $ leads to the following Dirac equation [53 



[7 M (i^-E^-(m-E.)]* = 
with the selfenergies S M and £ s given by 



(9) 



-v i d a +m 



— v % d a+m 

S s = g^ae a 



v i d a +m 



(10) 
(11) 



Both Lorentz-components of the selfenergy, S M and S g , show a linear dependence with 
respect to the meson fields p M , and a, as in the standard RMF. However, they contain 
an additional non-linear dependence on the partial derivatives. 

The following Proca and Klein-Gordon equations for the meson fields are obtained 



+ m\p v 

ri 2 



d u d"a + mla 



9w 

9a 



— i d a v a +m — -v a i d a+m 

\I/e a yip + \Ery e a \]> 



\]>e a 7^r\l> + $fy e A * 



— - idgv ±™ _ -=r- — '9n-f"i T 

\I/e ^ v[> + v]/e a v]/ 



(12) 
(13) 
(14) 



where Eqs. (fT2| [131 H4"|) show a similar form as in conventional RMF (A — > oo), except of 
the highly non-linear structure in the source terms. 
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3. Isospin-asymmetric nuclear matter 

We apply now the NLD model to isospin-asymmetric nuclear matter. In the RMF 
approximation all meson fields are treated as classical fields. In this approximation the 
spatial components of the vector meson fields vanish. Moreover, only the third component 
of the p-meson field in isospin space survives. To simplify the notation, the time-like 
component of the u° field will be labelled as u. Similar notation holds also for the time-like 
and third components in Lorentz- and isospin-space, respectively, of the vector-isovector 
p° field, which will be denoted as p. The baryon density will be labelled as ps- For the 
auxiliary vector we choose ir = (1,0). Note that, additional rearrangement terms can 
appear for different choices of tr in the NLD formalism. However, they exactly vanish in 



the mean-field approximation to infinite nuclear matter [53 



The Dirac equation, see Eq. (J5J), can be solved using the following ansatz 

*(s,p) = u(s,p)e- i '* x »( *) , (15) 



n 



where s and p = (E,p) stand, respectively, for the spin and 4- momentum of the proton (p) 
or neutron (n), denotes the space-time coordinate, and u is the Dirac spinor for positive 
energy states. The Dirac equation ([9]) maintains its form in infinite nuclear matter with 
selfenergies given by 



El = E vi = g w u e 'a ±g p pe a , (16) 

E^—m 

^si = 9*° e a— . (17) 



The upper (lower) sign in Eq. ( fl6l) refers from now on always to protons, i = p, (neutrons, 
i = n). The Dirac equation is solved inserting Eqs. (TT5]) : (TIB"]) and (TIT]) into Eq. flS]) for 
protons (i = p) or neutrons (i = n) 

rr E*% = (i-p + m*)%, (18) 

with the in-medium energy and mass given by 

E* :=Ei-( g u ue a ± g ppe ~ , (19) 



and 

m* := m — g a cre ~ . (20) 
These quantities are related to each other through the dispersion relation 

Ef-p 2 = mf , (21) 

which, because of the p-field, is different for protons and neutrons. The solutions of the 
Dirac equations for protons and neutrons are found in the usual way 

»,(>■/') - A', ( sj S , I . (22) 

E*+m* " : 
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with the spin eigenfunctions 4> s . The factor Ni is defined as Ni = ■\J E * 2 E*' i an d normalizes 
the spinors as follows Ui(s,p)^ Q Ui(s,p) = 1 and Ui(s,p)ui(s,p) = 

i 

In nuclear matter the meson field equations of motion read 

m l° = 9aps , m 2 p P = 9 P Pi , mlu = g^po , (23) 

where the vector-isoscalar p , vector-isovector pj and scalar-isoscalar p s source terms take 
the forms (k = 2 is the spin degeneracy factor) 

E i -m ft * I" JJl* _E i -m 

(^e-— ^) = — 2^ j d'p^e —&(p-p Ft ) = p sp + Psn , 

i=p,n i=p,n 1 

_ E^ — m r — ^ / Q E. i ~ — 

(^a e = T^p 2^ / d pe ~ e (P -PfJ = Pop + Pon , 

i=p,n i=p,n 

n E^—m 

(^a r 3 e — *») = Po P - Pon • (24) 

j=p,n 

Contrary to the Walecka model, in the NLD approach the vector-density po is not a con- 
served nucleon density. Latter one has to be derived from the generalized Noether-theorem, 
see Eq. Following Ref. j53j the conserved nucleon density after the resummation pro- 
cedure takes the form 

J° = Pb = £ [<^7°*i> (25) 

i=p,n 

(]„ Ei — rn O, . r. Ei—m C] „ Ea—tu 

- ^(^ e -^^)a + ^(^7° e -^^)a; + |(vl/ i7 e -^r3^)p , 

where the expectation value pw = J2i= P n(^7°^«) * s J us ^ the usua l density of the Walecka 
model [56 . 

The equation of state (EoS) is obtained from the 00-component of the energy-momentum 
tensor T^ u . The resummation procedure described in [53j results in 



i=p,n 



i=p,n i=p,n i=p,n 

+ l -{mla 2 -mlu 2 -my) , (26) 
with the additional expectation values 

Pf - 



i=p,n ^ ' i=p,n ' 

„ E^ — m K. w — -\ In Ea—ttl 

(%E i7 °e-^%) = _ £ / tPpEte-^eip - p Fi ) , 

i=p,n i=p,n 

Pf = E^TV^Ta^Hp* -p*. (27) 



The pressure P is given as usual by 



3 

However, since in the NLD approach rearrangement terms do not appear 53|, thermody- 
namic consistency is always fulfilled. Then, the pressure can be obtained from the energy 
density by using the thermodynamic relation 

P-Pl 9 -^ . (29) 

OPB 

Note that in Eq. ( 1291) the conserved baryon density ps from Eq. ( 1251) appears, and not the 
standard Walecka density. 

Two important effects appear in NLD model, which are absent in the minimal aup- 
models. At first, the selfenergies depend not only on density, but also explicitly on single- 



particle energy, as discussed in detail in Ref . 53] . These feature is retained also in isospin- 



asymmetric nuclear matter. As an additional novel feature of the NLD model, not only 
the Lorentz-vector component of the nucleon selfenergy, but also its Lorentz-scalar part 
depends implicitly (through the dispersion relation) on the isovector p-meson field. This 
feature arises from the isospin dependence of the single-particle energy, and it will generate 
an in-medium mass splitting between protons and neutrons in isospin asymmetric matter. 
Note that the a field is also affected by isospin effects, but the residual impact of the 
dispersion relation is negligibly small in this case. 

4. Numerical realization 

The meson-field equations depend explicitly on the baryon density, but only implicitly 
on the single-particle energy. Latter quantity is integrated out in the source terms. There- 
fore, it is sufficient to solve selfconsistently only the coupled set of meson-field equations 
(|23l) at a given density, and then determine the proton or neutron single-particle energy by 
solving numerically the dispersion relation (|2TT) for a given proton or neutron momentum 
relative to the nuclear matter at rest. 

For the density dependence we proceed as follows: the input parameters are the Fermi- 
momentum and the isospin-asymmetry. The isospin-asymmetry parameter a is determined 
according to the standard definition 

7° - 7° 

The set of the meson field equations ff23l) is solved then selfconsistently until convergence 
is achieved. For each iteration step the integrals in Eq. f f24|) are evaluated numerically, 
where at each momentum step the dispersion relation, see Eq. (I2T1) . has to be calculated 
also numerically. The whole procedure provides the meson-fields, with which we calculate 
the corresponding conserved baryon density according Eq. (125]) and the isospin-asymmetry 



8 



parameter according Eq. ([30]) . Note again, that the conserved baryon density ps differs 
from the Walecka one p w . Finally, the single-particle energies, which appear explicitly in 
the selfenergies, are evaluated at the corresponding proton or neutron Fermi-momentum. 

For the energy dependence at a fixed baryon density the procedure is similar as outlined 
above. We first calculate the meson-fields at a given density and isospin-asymmetry. Then 
we solve numerically the dispersion relation for protons (neutrons). In the first case the 
additional input parameter is the proton (neutron) momentum by keeping the momentum 
of the other isospin state at the corresponding fixed Fermi-momentum. The in-medium 
energy of a proton (i = p) or neutron (i = n) is defined then as 



Ei = ^Jmf+pZ + Ki , (31) 

where the proton or neutron Lorentz- vector selfenergy is given by Eq. (fTS]) . 

In principle, the NLD model contains no parameters. The original isoscalar aNN and 



uNN couplings can be taken from any linear Walecka model, e.g., [55[, as it has been 



done in Ref. [53| and in this work. The pNN coupling g p = 3 close to its universal value is 



used [60J which is also well determined by the one-boson-exchange model [581 ] . A small (big) 
numerical value of A provides a stronger (weaker) density and energy dependence of the 
mean-field, which influences the density dependence of the equation of state for densities 
at and above saturation and, at the same time, strongly affects the energy dependence 
of the Schrodinger-equivalent optical potential at saturation density. The simplest way 
followed here is to fix A just from the empirically well-known energy dependence of the 



Schrodinger equivalent optical potential. The cut-off parameter A is set to 770 MeV [53 
A different numerical value for the cut-off parameter is of course possible if one uses a 
different Walecka-like parametrization than the original one or if one adopts a different 
functional form for the non-linear operator than the simpler exponential one. 



5. Results and discussions 

Before presenting the calculations for isospin asymmetric nuclear matter we discuss 
briefly the NLD results for the special cases of symmetric matter (a = 0) and pure neutron 
matter (a = 1). In Fig. [I] we present the EoS for a = and a = 1 in terms of the 
pressure density as function of the baryon density. Here we compare the empirical heavy- 
ion data [ill ], see the shaded bands in Fig. [TJ with NLD and Walecka models. The inserted 
panels in Fig. [T] shows the nuclear EoS in terms of the binding energy per nucleon. The 
conventional Walecka model (dashed curves) leads to a very stiff EoS for symmetric nuclear 
matter and also pure neutron matter with a high compression modulus at saturation, 
because of the well-known linear divergent behaviour of the repulsive u field. On the other 
hand, the Lorentz-vector mean-field is suppressed to large extend at high baryon densities 
in the NLD model. See further discussions in Ref. [53[. This effect, which is retained 
also for the extreme case of pure neutron matter, weakens considerably the stiffness of the 
high-density EoS, see the insert in Fig. [TJ and the empirical region in the pressure-density 
phase diagram is reproduced fairly well by the NLD approach (solid curve). 
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baryon density p R [fm ] 



baryon density p R [fm ] 



Figure 1: The nuclear EoS for symmetric nuclear matter (main panel on the left) and pure neutron 
matter (main panel on the right) in terms of the pressure density in the Walecka (dashed curve) and 
NLD (solid curve) models. The shaded bands refer to the experimental regions extracted from heavy-ion 
collisions 0, The inserted panels on the left (right) describe the nuclear EoS for symmetric matter 
(left) and neutron matter (right) in terms of the binding energy per nucleon Eb — c/pb — m, again for 
both NLD and Walecka models. 



However, not only the density dependence of the nuclear mean-field, but also its energy 
dependence matters for nuclear reactions. This can be demonstrated using the in-medium 
optical potential which is defined as follows 

U opt = -E„ - S s + J- - T? v ) . (32) 

As shown in Fig. 121 the empirical data from Dirac phenomenology (open diamonds) pre- 
dict a saturation of the in-medium nucleon optical potential j42[ at intermediate energies. 
As well known, standard Walecka-type RMF models (RHD, dashed curve) are not able to 
reproduce this empirical saturation, and strongly overestimate the data at high energies. 
This disagreement remains also in other RMF models (DD, dot-dot-dashed curve), where 
the meson-nucleon couplings depend explicitly on the baryon density. A better description 
can be achieved here by the consideration of additional energy dependent terms in the in- 
teraction Lagrangian (D 3 C, dot-dashed curve). However, with the cost of many additional 
parameters [37|. Microscopic DBHF models (DBHF, dot-dashed-dashed curve) reproduce 



the optical potential in average up to energies close to pion production threshold [59|, [62 
On the other hand, the NLD model (NLD, solid curve) not only softens the EoS at high 
densities, but simultaneously generates the correct energy behaviour of the nucleon-nucleus 
optical potential at high energies. Note that, the same NLD approach describes the puz- 



zling energy dependence of the antiproton-nucleus potential, see Ref. [54 
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200 




kinetic energy E [MeV] 



Figure 2: Kinetic energy dependence Ekin = E — m of the Schrodinger equivalent optical potential at 
ground state density p sa t — 0.16 fm~ 3 in the linear Walecka (RHD, dashed curve) and NLD (solid curve) 
models. The other theoretical calculations, DD (dot-dot-dashed curve), D 3 C (dot-dashed curve) and 
DBHF (dot-dashed-dashed curve), are taken from Refs. [37|,|62|. All model calculations are compared with 
data from Dirac phenomenology (open diamonds) [421 ] . 

In Fig. [3] we present the density dependence of the nucleon selfenergies for various 
asymmetry parameters a. In symmetric nuclear matter (a = 0), the NLD vector 
and scalar S s selfenergies (solid curves) saturate at high densities. As one can see, they 
essentially follow the density dependence of the microscopic DBHF calculations [Hif (filled 
squares). As discussed in detail in [53], this saturation of the selfenergies with increasing 
density, in particular, of the repulsive vector part, is responsible for a softening of the EoS 
at high compressions. The saturation in density remains in the general case of isospin- 
asymmetric nuclear matter, see the dashed (a = 0.5) and dot-dashed (a = 1) curves in 
Fig. [3j As a novel feature, not only the vector, but also the scalar component of the nucleon 
selfenergy gets isospin dependent. The isospin dependence of the scalar selfenergy, observed 
here, arises from the p-field induced isospin splitting between the in-medium proton and 
neutron single particle energies. This effect will be discussed in more details later on. 

The saturation of the selfenergies in density affects considerably the density dependence 
of the symmetry energy. The symmetry energy is extracted from the empirical parabolic 
low and defined as the second derivative of the energy per nucleon versus the asymmetry 
parameter a, or as the energy difference between pure neutron matter and symmetric 
nuclear matter. Fig. H] (left panel) shows the density dependence of the symmetry energy 
within the NLD (solid curve) and the Walecka (dashed curve) models. The symmetry 
energy E sym = T sym + V sym in the standard RMF grows almost linearly with baryon 
density. Here the symmetry potential V sym is proportional to the p field, which increases 
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Figure 3: Density dependence of the vector £„ (upper panel) and scalar E s (lower panel) selfenergies. The 
NLD results for symmetric nuclear matter (solid curve) arc compared with DBHF calculations 63] (filled 
squares). The NLD fields for a — 0.5 (dashed curves) and a = 1 (dot-dashed curves) are shown too, and 
the direction of the proton and neutron splitting is marked by the arrows, as indicated. 



linearly with density, pg. The moderate saturation at very high densities arises from the 
suppression of the kinetic part T sym ~ ^W- (with E* F being the in-medium Fermi-energy). 

F 

The NLD model gives a different behaviour at high baryon densities with respect to the 
conventional RMF. A considerable softening of the high density EoS is obtained, but now 
for the isovector sector. 



Available data for the symmetry energy exist around saturation density only [64( . The 



theoretical models describe the experimental data fairly well, however, they give different 
predictions for the high-density symmetry energy, see for reviews Refs. |2j. In non- 
relativistic models the adjustment of the model parameters to the empirical saturation 
point does not necessarily constraint the stiffness of the symmetry energy. In fact, a soft and 
a stiff symmetry energy can be obtained at high densities [2J. In the relativistic Walecka- 



type models, the high density dependence of the symmetry energy is always stiff [3l|, [32 



Similar results are obtained in the DBHF approach and within density-dependent rela- 
tivistic Hartree-Fock (RHF) approaches 40|, |4l| , see the dashed and dot-dot-dashed curves 
in the right panel of FigJU where the symmetry energy becomes stiffer with increasing 
density [38|. It is interesting to note that the RHF models contains also an energy de- 
pendence, apart the density dependent coupling functions. On the other hand, the DD 



models |37|, |36| , see the dot-dashed (DD) and dot (DD-ME5) curves in the right panel 
of FigJU predict a rather soft symmetry energy, which is similar to the NLD results (solid 
curve) . 

Direct experimental access to the high-density region of the symmetry energy is still 
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Figure 4: (Left panel) Density dependence of the symmetry energy -Esym in the standard linear Walecka 
model (dashed curve) and NLD approach (solid curve). (Right panel) The same as in the left panel, but 
in comparison between the NLD model with other theoretical approaches: DBHF (dashed curve) 13811 . DD 
(dot-dashed curve) [HEl, RHF (dot-dot-dashed curve) [H| and DD-MEJ (thick-doted curve) (Sf. The 
symbols refer to experimental data from Rcf. 64]. 



incompletely known. A plausible way to constraint the stiffness of the symmetry energy 
at supra- normal densities consists in studies of compact neutron stars [4]. In fact, the 
symmetry energy must saturate or, at least, show a rather soft behaviour at high baryon 
densities in order to respect the direct URCA-limit [4j]. This constraint is satisfied in the 
NLD model. 



The density dependence of the nucleon effective or Dirac mass in the NLD model [53 



is in remarkable agreement with the DBHF calculations. The insert in Fig. [5] describes the 
Dirac mass of the nucleon in symmetric nuclear matter. Here the solid curve corres pon ds 



to the NLD results and the square symbols describe the DBHF calculations of Ref. [63 . 

Interestingly, in the NLD model an effective mass splitting between protons and neu- 
trons is generated in asymmetric matter by taking into account at the Hartree-level the 
p-field only. This is demonstrated in Fig. [5]for two asymmetry parameters a = 0.5 (dashed 
curves) and a = 1 (dot-dashed curves). In particular, a mass splitting m* > m* is ob- 
tained for neutron-rich matter, which becomes pronounced with increasing density and 
isospin asymmetry parameter. The in-medium mass splitting in the NLD model is due 
to the isospin dependence of the single particle energy which appears in the selfenergies. 
Note that, any explicit energy dependence of the scalar selfenergy will result in a p- induced 
mass splitting of nucleons in isospin asymmetric nuclear matter. 

As discussed in the introduction, an effective mass splitting between protons and neu- 
trons do not occur in the Walecka-type cxcup-models solved at the Hartree-level. The 
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Figure 5: Density dependence of the Dirac mass splitting between protons and neutrons within the NLD 
model. The asymmetry parameters are a — (solid curve), a = 0.5 (dashed curves) and a — 1 (dot-dashed 
curve). The curves above (below) the solid curve (symmetric matter) correspond to the neutron (proton) 
effective masses, as indicated. The insert describes the Dirac mass of the nucleon in symmetric nuclear 
matter in the NLD (solid) and DBHF (filled squares) calculations. 



introduction of the scalar-isovector 5-meson in the nuclear mean-field is necessary for this 
However, one gets more parameters to fit the bulk asymmetry parameter 



31, 32 



purpose 

at saturation density and, in particular, one needs to re- adjust again the pNN coupling. 
This procedure can be ambiguous, since the SNN coupling is not fixed neither by nuclear 
structure 34[ nor by the one-boson-exchange models 



63 . On the other hand, as it has 



been investigated intensively in the past, the 5- meson implies interesting effects if one goes 
beyond the saturation density of ground state nuclear matter and study its influence in 



compressed matter created in heavy-ion collisions [24l. |25 



perimental data do not allow still decisive conclusio ns [65| 



Unfortunately, the available ex- 
The cxa;p<5-models predict a mass 
splitting m* < m* for neutron-rich matter [l], EH 32, 24 1. However, just from the energy 
dependence of the scalar selfenergy the NLD model predicts an opposite splitting m* > m* 
for neutron-rich matter. This trend is unique and cannot be reversed in the NLD model. 
We remind here, that RMF models treated on the Hartree-Fock level (RHF) 0, 41 



in which the selfenergies depend also explicitly on the single-particle energy due to the 
Fock-terms, predict for the non-relativistic effective mass a mass splitting m, 
densities p# < 0.8p sat , but this trend is reversed for densities p# > O.Sp sat 



40 



> m* up to 
Therefore, 

the property m* < m* is not a genuine feature of the RMF models, as has been believed 
so far. 

The energy dependence of the isovector nuclear potential is described by the Schrodinger- 
equivalent optical potential, see Eq. (I3"2"j) . but now separated between protons and neutrons. 
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Figure 6: Momentum dependence of the Schrddinger equivalent neutron (left panel) and proton (right 
panel) optical potentials at various asymmetries, as indicated. Calculations in the NLD model are shown 
only. The standard RHD model does not contain a momentum dependence. 



Fig. E] shows these potentials as function of the corresponding particle momentum and 
different asymmetry parameters a. The observed isospin dependence of the proton (right 
panel) and neutron (left panel) optical potentials results from cancellation effects in Eq. ( 132]) 
between different isospin dependent components of the selfenergies. It is remarkable that 



the microscopic (D)BHF calculations [38|, |48|, |49[ predict very similar behaviour as the NLD 



model. This is shown in Fig. where the proton and neutron optical potentials are shown 
as a function of the asymmetry parameter a at fixed density p = p sat and three different 
momenta, as indicated in Fig. [7J Indeed, at low momenta up to the Fermi momentum (left 
and middle panels) the isospin splitting between the proton and neutron optical potentials 
in the NLD model (solid and dashed curves) follows closely the microscopic calculations 
(filled and open squares for DBHF, filled and open triangles for BHF). At higher momenta 
(right panel) the NLD isospin splitting is similar to the (D)BHF calculations, up to an 
absolute value. This is related to a stronger energy dependence in the NLD model relative 
to these microscopic approaches. Recall, that this strong energy dependence of the NLD 
potential is essential for the description of the Schrodinger-equivalent optical potential 
shown in Fig. [2J 

The energy dependence of the isovector optical potential is conventionally described by 
the Lane-type potential. It is defined as follows [H, 0] 

t4o = , (33) 

where U p ^ n denotes the proton/neutron Schrodinger-equivalent optical potentials. Fig. E] 
shows the momentum dependence of the Lane-potential in the NLD approach (solid curve). 
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Figure 7: Asymmetry parameter a dependence of the optical potentials at fixed saturation density and at 
three fixed momenta (p = 0, 1, 2 fm^ 1 as indicated) in the NLD (solid and dashed curves for neutrons and 
protons), in the DBHF [38| (filled and open squares for protons and neutrons) and in the BHF 48] (filled 
and open triangles for protons and neutrons) models. The arrows in the left panel indicate the direction 
of the proton and neutron splitting in these three figures. 



Also the results of linear Walecka model denoted as RHD (dashed curve) and different 
microscopic approaches are shown. The linear Walecka model predicts an almost linear 
behaviour in energy, as in the isoscalar case, and it diverges for p — > oo. The Lane potential 
in the NLD model shows an opposite behaviour with respect to RHD. It decreases with 
increasing particle momentum. In contrast to the symmetric proton- and antiproton- 
nucleus optical potentials 53|, |54(, the empirical situation in the isovector case is still 
not well known. Different empirical studies predict an opposite behaviour of Ui so versus 
momentum. In fact, the analysis of Madland et al. 4fJ, |47| predicts an increase of the 



isovector potential with increasing momentum, while in the analysis of Lane et al. [43|, |44 



45 



a decrease of U\ so with rising momentum has been obtained. Nevertheless, as shown in 
Fig. [HI the NLD calculations agree with the calculations from the microscopic Brueckner 
theory which also show the decrease of the Lane potential as a function of the momentum. 



6. Summary 

In summary, the NLD formalism which incorporates on a mean-field level the energy 
dependence of the nucleon selfenergies has been applied to isospin-asymmetric nuclear 
matter. The model describes the bulk properties of symmetric matter and compares well 
with the results from microscopic Brueckner calculations. Due to the explicit energy de- 
pendence of the nuclear mean-field, the NLD model reproduces very well the empirical 
proton and antiproton optical potentials [53|, |54|. Same energy dependence is responsible 
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Figure 8: Energy dependence of the Lane potential at saturation density and an isospin asymmetry 
parameter a — 0.4. The NLD model (solid curve) is compared with the Walecka model (RHD, dashed 
curve) and with the DBHF [3|| (solid curve with filled squares) and BHF [49[ (dot-dashed curve with 
filled triangles) models. 

for the softening of the EoS at high densities. In the present model the isospin effects 
are generated by using the conventional p-meson exchange. Such a minimal extension of 
the NLD formalism describes the symmetry energy around saturation and predict its soft- 
ening (saturation) at high baryon densities. This feature is an important element in the 
description of dense nuclear systems created in heavy ion collisions and/or in the interior 
of neutron stars. However, such a soft symmetry energy is at variance with standard RHD 
models, where it rises linearly with increasing density. 

As a novel feature we find that the energy dependent NLD interactions generate a mass 
splitting between protons and neutrons already within the conventional p-exchange scheme. 
In particular, a Dirac mass splitting m* > m* is obtained for neutron-rich matter. However, 
we point out that the p-induced nucleon mass splitting is a rather general feature of energy 
dependent interactions and should show up in any model which incorporates them. 

We further discussed the energy dependence of the isovector optical potentials. We 
found that the NLD model leads to a decreasing Lane-type potential with increasing parti- 
cle momentum. This disagrees with the conventional RMF models, but is consistent with 
the results from microscopic Brueckner calculations. Also the comparison of the isovec- 
tor splitting in the proton/neutron optical potentials agrees well with the results from 
DBHF and BHF microscopic approaches. These suggest that the NLD croup mean-field 
model which accommodates the energy dependence of nucleon selfenergies, is in qualita- 
tive agreement with microscopic (Dirac) Brueckner calculations. Since the NLD model 
describes also asymmetric matter features fairly well and agrees with microscopic DBHF 



17 



models, it could be interesting to apply the NLD approach to finite nuclei. However, due 
to the energy dependence of the fields, the main advantage of the NLD approach would be 
its application to the description of heavy-ion reactions within the transport theory. 
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